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ABSTRACT 


EV^r  two  locations  vri.thin  the  surf  zone  sea  surface  elevations  vrere 
daserved  using  a wave  staff  and  a pressure  sensor  while  simultaneously  the 
•two  harizontal  orthogonal  ccmponentS/  u and  v»  of  water  particle  velocity 
were  measured. 

Surface  elevations  derived  fran  pressure  sensors  are  lower,  mainly  on 
the  region  of  the  crest,  ccnpared  vjith  the  same  surface  elevations  measured 
with  wave  gages.  Pressure  records  are  more  smoothed  than  wave  gage  records 
and  the  energy  ocnpited  for  waves  measured  with  a pressure  sensor  is  con- 
sistently snaller  than  for  waves  measured  with  a wave  gage. 

Methods  for  ccnverting  pressure  to  surface  elevation  are  given  which 
include  the  non-linear  velocity  term  (u^+v^)  which  is  usually  neglected  in 
the  Bernoulli  equation.  Two  techniques  are  proposed  to  include  this  term: 
1)  flowmeters  are  used  to  measure  u and  v,  and  2)  the  Bernoulli  term  is 
derived  by  determining  the  velocities  by  convolving  the  pressure  records 
using  a weighting  function  determined  frcm  shallow  vaster  theory. 

The  first  tecdmique  gave  ircsrovements  on  the  order  of  3.4%  to  7.8%  of 
the  total  variance  of  the  wave  gage  spectrum,  whereas  the  second  technique 
gave  improvements  on  the  order  of  3.1%  to  9.7%.  The  iitprovements  in  both 
^•^ag^ac  are  approximately  the  same,  with  the  second  technique  having  the 
advantage  of  requiring  only  a pressure  transducer. 
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I . INTRODUCTION 


A.  HISTORICAL  PERSPECTIVE 

In  studying  shallow  water  and  breaking  waves,  many  theo- 
retical and  practical  problems  are  encountered.  Analytical 
wave  theories  do  not  correctly  characterize  the  physics  of 
wave  breaking.  The  environment  is  neither  friendly  nor  easily 
reproduceable  in  a laboratory.  A primary  problem  that  is  en- 
countered is  how  to  measure  the  wave  height  at  or  near  the 
breaking  point.  In  general,  the  two  basic  types  of  sensors 
used  to  measure  the  surface  elevations  have  been  subsurface 
pressure  sensors  and  surface  piercing  wave  gages. 

The  use  of  pressure  sensors  to  infer  surface  elevation  is 
particularly  desirable  as  a means  of  conveniently  measuring 
breaking  waves,  especially  high  energy  waves.  Pressure  sensors 
have  been  used  more  commonly  because,  compared  with  surface 
piercing  gages,  they  are  rugged,  less  vulnerable  to  environ- 
mental hazard,  and  easier  to  install;  in  many  locations,  pres- 
sure sensors  are  the  only  feasible  means  of  measuring  the  waves 

The  conventional  procedure  used  to  infer  surface  elevation 
from  a pressure  signal  is  to  calculate  the  energy  density  spec- 
trum of  a pressure  record  and  to  apply  a spectral  transfer  func 
tion  derived  from  linear  theory.  The  pressure  pulse  is  atten- 
uated with  depth,  with  the  attenuation  increasing  with  wave 
frequency.  Hence,  the  spectral  transfer  function  is  dependent 
on  frequency  and  depth. 


Esteva  and  Harris  (1970)  compared  results  from  a wave  staff 
and  a pressure  transducer  in  16  feet  of  water.  The  pressure 
spectrum  was  compensated  using  linear  wave  theory  to  obtain  the 
transfer  function.  The  computed  surface  energy  density  spectrum 
was  then  summed  over  all  components  and  the  compensated  root 
mean  square  of  the  total  energy  determined;  this  procedure  led 
to  good  estimates  of  wave  heights  from  pressure  records  as  com- 
pared with  those  determined  from  surface  records.  However,  a 
number  of  investigators  have  found  discrepancies  using  the 
linear  transfer  function  in  intermediate  to  shallow  water. 
Takahashi  et  al.  (1967)  , Glukhovski  (1961)  , Shooter  and  Ellis 
(1967) , Gerhardt  et  al.  (1955) , suggest  several  different  multi- 
plicative correction  factors  to  the  linear  transfer  function 
for  various  conditions.  In  general  the  correction  factors  de- 
crease with  decreasing  period,  being  greater  than  1.0  for  long- 
period  waves  and  less  than  1.0  for  short-period  waves  (Shore 
Protection  Manual,  U.  S.  Army,  Corps  of  Engineers,  1975). 

Near  or  at  breaking,  linear  wave  theory  no  longer  is  appli- 
cable, and  using  the  linear  transfer  function  does  not  give  as 
good  an  approximation.  Van  Dorn  (1977) , in  a laboratory  study 
of  breaking  waves,  found  that  surface  elevations  computed  from 
pressure  sensors  were  correctly  measured  in  the  troughs  but 
were  substantially  less  under  the  wave  crests.  At  the  crests, 
the  surface  elevations  were  sometimes  50  per  cent  too  low  as 
compared  with  wave  staff  measurements  (Figure  1) . The  discre- 
pancy between  pressure  and  wave  staff  measurements  is  caused 
I 

t by  a dynamic  pressure  reduction  due  to  the  increased  velocities 
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ELEVATION,  cm 


DEPTH,  D-cm 


Figure  1.  Hodographs  of  maxiaium  and  minimum  waves 
elevations  as  an  individual  wave  moves 
toward  shore.  Pressure  and  staff  elevations 
agree  under  troughs,  o-t  F.«a<  pressure  is 
clearly  a poor  indira*.-  --f  -rest  elevation 
(S  is  beach  slope 
From  Van  Dorn  1^' 
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under  the  crests  of  breaking  waves.  Hence,  the  pressure 
records  are  much  smoother  and  rounded  off  at  the  crests  as 
compared  with  the  surface  staff  measurements  (Thornton  et  al., 
1976) . 

B . OBJECTIVE 

The  objective  of  thi,'?  research  is  to  improve  techniques 
for  inferring  surface  elevations  from  the  pressure  sensors 
records  and  to  test  the  techniques  in  the  field  under  break- 
ing waves. 

Pressure  is  converted  to  surface  elevation  by  including 
the  non-linear  velocity  term  which  is  usually  neglected  in  the 
Bernoulli  equation.  Two  techniques  are  proposed  for  including 
the  Bernoulli  term: 

1.  The  Bernoulli  term  can  be  measured  directly  using 
flowmeters . 

I 2.  The  Bernoulli  term  can  be  derived  by  convolving  the 

pressure  records  to  determine  the  velocities  using 
shallow  water  linear  theory.  This  technique  has  the 
advantage  of  requiring  only  pressure  sensors. 


Both  techniques  are  equivalent  to  a second  order  perturbation 
scheme  v/hich  includes  the  neglected  non-linear  term. 


II.  MEASUREMENTS 


f 

A.  EXPERIMENT  SITES 

For  two  locations  within  the  surf  zone  sea  surface  eleva- 
tions were  observed  using  a wave  staff  and  a pressure  sensor 
while  simultaneously  the  two  horizontal  orthogonal  components, 
u and  V,  of  water  particle  velocity  were  measured.  The  data 
used  in  this  paper  were  measured  at  Del  Monte  Beach  and  Scripps 
Beach,  California. 

Experiments  were  conducted  at  Del  Monte  Beach  within 
Monterey  Bay  on  8 March  1978.  The  waves  at  this  location  were 
plunging-spilling  breakers.  Because  of  topographic  sheltering 
and  severe  directional  filtering  due  to  refraction  by  the 
geometry  of  the  bay,  the  waves  offshore  were  narrow  band  swell 
type  waves  which  impinged  perpendicular  to  the  shore;  hence, 
a simplification  to  a two  dimensional  narrowbanded  wave  de- 
scription is  allowed.  The  median  grain  size  is  approximately 
0.2  mm  (taken  at  the  watar  line)  and  the  beach  slope  varied 
between  1:14  and  1:40. 

Experiments  were  conducted  at  Scripps  Beach  at  La  Jolla 
on  19-20  July  1978,  during  the  morning  and  evening  of  each  day. 
Scripps  Beach  has  an  approximately  planar  1:80  slope.  The 
median  grain  size  was  0.1  mm.  Spilling  breakers  predominated. 


A typical  beach  profile  for  each  experiment  is  shown  in 
Appendix  A-. 


B . INSTRUMENTATION 


During  the  Del  Monte  Beach  experiments,  the  instruments 
used  were  a capacitance  wave  gage,  two  Marsh-McBirney  model 
721  electromagnetic  current  meters  and  a Statham  model  PA506 
pressure  transducer. 

Measurements  at  Scripps  Beach  were  made  using  a resistance 
wire  wave  gage,  one  Marsh-McBirney  model  512  flowmeter  and  the 
same  pressure  sensor  used  at  Del  Monte  Beach. 

1.  Capacitance  Wave  Gage 

The  capacitance  wave  gage  was  fashioned  from  3/8-inch 
outside  diameter  stainless  steel  rod.  The  rod  was  tightly 
covered  with  1/16-inch  wall  thickness  polypropylene  tubing. 

The  gage  operates  on  the  principle  that  a change  in  the  plate 
dimension  of  a capacitor  changes  its  capacitance  and  conse- 
quently an  output  voltage.  The  insulated  steel  rod  and  sea- 
water act  as  the  plates,  the  insulation  functioning  as  the 
dielectric.  As  the  surface  elevation  fluctuates,  the  capaci- 
tance of  the  circuit  changes,  and  the  output  voltage  responds 
linearly.  These  output  voltage  fluctuations  were  sensed  by  a 
transistorized  circuit  powered  from  the  beach.  The  circuit 
was  designed  by  McGoldrick  (1969) . The  electronics  package 
was  housed  in  a watertight  brass  case  which  was  mounted  on  a 
tower  during  the  8 March  experiment.  This  allowed  the  con- 
necting leads  to  be  less  than  30  cm  long,  thereby  minimizing 
wire-to-wire  capacitance.  The  gage  was  statically  calibrated 
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in  the  laboratory  prior  to  the  experiment.  Accuracy  was 
estimated  to  .005  m.  For  all  experiments , the  calibration 
factors  are  shown  in  Appendix  B. 


2.  Resistance  Wave  Gage 

The  resistance  wave  gage  is  a surface-piercing,  dual- 
wire, resistance  sensor.  The  wave  staff  wires  were  mounted  on 
a 1.5-inch  outside  diameter,  1/8-inch  wall  fiberglass  fishing 
pole  of  5-m  length,  epoxied  into  an  aluminum  flange  at  the 
bottom.  The  two  resistance  wires  are  22  gage  nichrome,  attached 
to  a PVC  collar  at  the  top  of  the  pole  and  to  the  flange  at  the 
bottom.  The  wires  are  4 inches  apart  over  their  entire  length. 

A twisted-pair  cable  running  inside  the  fiberglass  pole  connects 
the  upper  ends  of  the  wires  and  provides  the  electrical  con- 
nection between  the  instrument  package  and  the  sensor  wires. 

The  gage  operates  on  the  principle  that  the  total  resistance 
measured  between  the  nichrome  wires  at  the  top  is  proportional 
to  the  length  of  the  exposed  (not  immersed)  wire.  The  wires 
are  the  unknown  resistance  in  an  AC  circuit,  whose  output  is 
linearly  related  to  the  length  of  the  exposed  wire. 

The  wave  gage  was  designed  and  built  by  the  Shore 
Processes  Laboratory  at  Scripps  Institution  of  Oceanography  at 
La  Jolla,  California;  field  and  laboratory  tests  have  shown  it 
I to  be  durable  both  inside  and  outside  of  the  surf  zone,  easily 

I installed  and  operated,  and  highly  linear.  The  instrument  has 

a resolution  of  a few  millimeters,  accuracy  better  than  one 
centimeter,  negligible  long-term  drift  and  excellent  temperature 
stability  and  frequence  response. 


3 . Flowmeters 

The  flowmeters  used  were  Marsh-McBirney  Models  721  and 
512  Electromagnetic  current  meters.  The  flowmeter  operation 
is  based  on  Faraday's  principle  of  electromagnetic  induction. 

Each  probe  measures  water  particle  velocity  in  two  orthogonal 
directions  through  a range  of  zero  to  three  meters  per  second,  ^ 

I 

with  a maximum  output  error  of  two  percent.  The  flowmeters 
were  dynaunically  calibrated  with  an  oscillating  platform  at- 
tached to  an  eccentric  aann  drive  by  a variable  speed  motor. 

Measurement  accuracy  was  determined  to  be  ±.02  m/sec  during 
calibration. 

t 

4.  Pressure  Transducer  I 

The  pressure  sensor  used  was  a Gould  Statham  Model 
PA506  thin  film  strain  gage.  The  sensing  element  is  a vacuum- 
deposited  resistive  balanced  fully  active  strain  gage  bridge 
with  an  output  voltage  responding  linearly  to  the  water  depth. 

The  sensor  measures  absolute  pressure  rather  than  relative, 
therefore  variations  in  atmospheric  pressure  occuring  during 
the  experiments  and  or  during  calibrations  will  slightly 
influence  the  apparent  mean  water  depth  measured.  The  instru- 

I 

i 

ment  has  a resolution  of  few  millimeters  and  a pressure  range 

from  3 to  20  psi.  The  pressure  sensor  was  statically  calibrated 

in  the  laboratory  prior  to  the  experiments  and  measurement  t 

accuracy  was  determined  to  be  ±.005  m.  ! 
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III.  ANALYSIS  OF  DATA 


Continuous  time  series  records  of  all  measurements  were 
collected  within  an  hour  or  two  on  both  the  ebb  and  flow  sides 
of  high  tide.  Analogue  records  were  made  on  both  magnetic 
tape  and  strip  charts. 

The  data  were  digitized  at  .2-sec  intervals  resulting  in 
a Nyquist  frequency  of  2.5  Hz  for  the  March  experiment  and  at 
.25-sec  intervals  resulting  in  a Nyquist  frequency  of  2.0  Hz 
for  the  July  data.  This  sampling  rate  was  sufficiently  high 
to  avoid  aliasing  of  energy  into  the  portion  of  the  spectra 
which  was  of  interest. 

A mean  value  was  computed  for  all  data  sets  and  the  data 
were  linearly  detrended  to  remove  tidal  effects.  Variance, 
standard  deviation,  skewness  and  kurtosis  of  the  distributions 
were  computed.  Because  the  record  lengths  were  slightly  dif- 
ferent for  each  run,  there  is  a small  difference  in  the  fre- 
quency resolution  for  each  run. 

The  power  spectra  of  the  wave  gage  and  pressure  sensor  were 
calculated  by  using  the  Fast  Fourier  Transform  (FFT)  technique. 
A section  of  the  record  approximately  35  minutes  long  was  used 
to  calculate  the  power  spectra.  The  35  minutes  were  chopped 
into  20  blocks.  Spectral  estimates  were  obtained  using  the 
FFT  were  ensemble  averaged  over  the  20  intervals  resulting 
in  approximately  40  degrees  of  freedom  for  each  spectral 


estimate. 


I 

j A data  window  was  applied  to  each  interval  to  minimize 

i 

leakage.  The  data  window  used  was  a cosine  square  taper  func- 
tion  over  the  first  and  last  5%  of  the  record. 

! Cross-spectra  were  computed  between  the  wave  and  pressure 

records  in  the  form  of  co-  and  quad-spectra.  Coherence  and 
phase  spectra  were  then  determined  from  the  cross-spectra. 
Coherence  and  phase  spectra  as  a function  of  frequency  indi- 
cate the  regions  and  degree  of  linear  relationship  and  phase 
between  the  wave  gage  and  pressure  sensor  data. 
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IV.  THEORY 


Assuming  an  irrotational , incompressible  flow,  a velocity 
potential  <{)  exists  and  is  described  as: 


(1) 


where  u,  v and  w are  the  velocity  components  in  the  horizontal 
and  vertical  directions  respectively.  Since  the  mean  hydro- 
static pressure  is  constant,  only  the  pressure  fluctuations 
about  the  mean,  Ap,  are  of  concern.  The  pressure  fluctuations 
are  described  by  the  Bernoulli  theorem  as: 

Ap  (t)  = p |1  - i p (u^  + v^  + w^)  (2) 

where  p is  density.  The  first  term  on  the  right  is  the  pres- 
sure due  to  the  wave  form  expressed  in  terms  of  the  velocity 
potential  and  the  second  term  is  the  contribution  to  Ap  due  to 
the  kinetic  energy  of  the  water  motion  induced  by  the  passage 
of  the  wave.  The  second  term  on  the  right  will  be  referred  to 
hereafter  as  the  Bernoulli  term. 

In  deep  water,  the  non-linear  Bernoulli  term  can  generally 
be  neglected  and  the  solution  for  linear  theory  is: 

Ap  (a,t)  = pg  ^ ri(CT,t)  = ri(a,t)  (3) 

cosh  kd  H(a) 
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= r ^ 

where  H(a)  is  the  linear  spectral  transfer  function  relating 
pressure  to  the  surface  n(t),  and  a is  the  radial  frequency. 

In  the  spectral  domain 

S^(a)  = |H(a)|^ 

In  shallow  water  and  in  the  vicinity  of  the  breaker  point, 

the  contribution  to  Ap  by  the  Bernoulli  term  becomes  important. 

Confining  the  discussion  to  shallow  water  waves,  w is  at  least 

2 

one  order  of  magnitude  less  than  u and  v (making  w very 

2 2 

small  compared  to  u and  v ) , so  that  the  contributxon  that  w 
makes  to  the  Bernoulli  term  can  be  neglected.  In  shallow  water, 
the  ratio  of  the  means  of  the  two  terms  on  the  right  hand  side 
of  (2)  with  pressure  meter  located  on  the  bottom  is  (Kinsman, 

1965)  : 

u^  + v^  2 (5) 

-TT- 

where  a is  the  amplitude  and  d is  the  depth.  Substituting 

the  saturation  relationship  at  breaking  for  solitary  wave  •'  | 

theory : . j 

a = .39  d (6) 

I 

into  (5)  shows  the  ratio  at  breaking  can  become  as  large  as  30%. 

Two  techniques  for  converting  pressure  to  surface  elevation 
are  proposed  which  include  the  Bernoulli  term.  Starting 

with  (2)  and  including  the  Bernoulli  term,  i 
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Ap  (a,t)  = — i—  n(a,t)  - (u^  + v^)  (7) 

H(a)  2 

or 

n(a,t)  = H(o)  [Ap(a,t)  + — (u^  + v^)  ] 

2 

= H(a)  f (a,t)  (8) 

In  the  first  technique  proposed,  u and  v are  measured  in 
the  field  with  electromagnetic  flowmeters  and  the  Bernoulli 
term  computed  and  added  to  the  differential  pressure.  The 
spectrum  of  the  surface  elevation  is  calculated  from  the 
combined  time  series 

(a)  = |H(a)  I ^ 3^(0)  (9) 

The  second  technique  requires  only  a pressure  measurement. 
The  velocity  used  in  the  Bernoulli  term  is  calculated  by  con- 
volving the  pressure  record 

u(t)  = h(T)  Ap  (t  - t)  (10) 

where  the  weighting  function  h(T)  is  given  by  the  Fourier 
transform  of  the  transfer  function 

= iT  ei2TrfT 

The  transfer  function  derived  from  linear  theory  can  be  used 
as  a first  approximation.  For  shallow  water  waves,  the  pres- 
sure is  essentially  hydrostatic  and  the  waves  are  non- 
dispersive.  Hence,  the  transfer  function  is  independent  of 
frequency  and  the  weighting  function  is  independent  of  time. 
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The  velocity  is  easily  calculated  for  this  simplified  case 


Ap  = c 


where  celerity,  c , is  given  as  a first  approximation  by: 


V.  RESULTS 


A.  QUALITATIVE  AND  QUANTITATIVE  DESCRIPTION  OF  THE 

OBSERVED  DATA 

Due  to  the  similarity  of  the  spectral  shapes  of  the  various 
analysed  data  the  discussion  of  results  will  be  exemplified  by 
the  data  of  19  July  evening.  Except  as  w.ere  noted,  the  19 
July  evening  results  are  representative  of  the  other  data. 
Appendix  C contains  the  spectra  calculated  for  all  experiments. 

Surface  elevations  measured  using  the  wave  staff  are  com- 
pared in  Fig.  2 to  the  hydrostatic  pressure  head  and  the 
hydrostatic  pressure  head  including  the  Bernoulli  term. 

The  pressure  sensor  and  wave  staff  records  generally  agree 
in  areas  near  the  wave  through,  but  the  pressure  sensor  under- 
estimates the  region  near  the  crest  as  shown  in  Fig.  2.  The 
results  are  in  accord  with  the  laboratory  study  of  Van  Dorn 
(1977).  Referring  to  Eq.  2,  the  hydrostatic  pressure  fluctua- 
tions under  crests  and  troughs  is  given  by: 

Ap  = pgn  “ ^ P (14) 

The  second  term  on  the  right  hand  side  of  Eq.  14,  i.e.,  the 
Bernoulli  term,  is  always  a reduction  in  pressure  relative 
to  hydrostatic;  since  u and  v become  maximum  under  crests  and 
minimum  at  the  troughs,  the  Ap  term  is  more  reduced  under  the 
crests  than  at  the  troughs  relative  to  hydrostatic. 


The  sea  surface  given  by  the  wave  staff  and  the  hydrostatic 
pressure  head  were  found  to  have  positive  skewness  and  kurtosis, 


with  larger  values  for  the  wave  staff  distribution.  Larger 
values  of  positive  skewness  imply  that  the  wave  staff  indicates 
larger  heights  for  the  crests.  Larger  kurtosis  for  the  wave 
staff  indicates,  as  expected,  more  peaked  crests  than  observed 
by  the  hydrostatic  pressure  head. 

The  variance,  standard  deviation,  skewness  and  kurtosis 
are  summarized  in  Table  I. 

A measure  of  the  total  energy  in  the  measured  spectra  is 

given  by  the  variance.  Variance  was  computed  for  all  data 

2 

sets.  For  the  wave  staff  it  ranged  from  .0335  to  .0471  m , 

2 

and  for  the  pressure  sensor  from  .0247  to  0.377  m . The 
variances  of  the  wave  staff  were  always  larger  than  for  the 
pressure  transducer. 

A mean  water  surface  elevation  was  computed  for  all  data 
sets.  The  location  of  the  breaker  point  is  highly  dependent 
on  the  depth  of  the  water.  The  tidal  range  at  both  locations 
is  approximately  1.5  meters,  resulting  in  the  breaker  location 
changing  considerably  during  a tidal  cycle.  During  a tidal 
cycle  it  was  possible  to  measure  waves  both  inside  and  outside 
the  surf  zone.  The  mean  depths  during  measurements,  calculated 
using  wave  staff  and  pressure  sensor  are  given  in  Table  II 
along  with  the  relative  location  of  the  mean  breaker  line  and 
type  of  breaker.  The  mean  sea  surface  values  measured  with 
the  wave  staff  ranged  from  1.57  to  2.03  m,  whereas  the  pressure 
sensor  measurements  ranged  from  1.39  to  1.87  m.  For  all  data. 
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TABLE  I.  VARIANCE,  STANDARD  DEVIATION,  SKEWNESS  AND  KURTOSIS  FOR  ALL  DATA 


T 


The  mean  water  depth  given  by  the  pressure  sensor  is  less  than 
the  mean  water  depth  measured  with  the  wave  staff  by  a differ- 
ence that  ranged  from  .08  to  .20  m.  (Figure  3 the  mean  water 
level  measured  using  wave  staff  and  pressure  sensor.)  This 
difference  can  be  explained  due  to  two  main  factors: 

1.  The  hydrostatic  pressure  head  is  decreased  by  the 
Bernoulli  term  and  therefore  gives  a mean  value  less 
than  the  mean  value  measured  by  the  wave  staff. 

2.  The  pressure  transducer  used  measures  absolute  pres- 
sure rather  than  relative;  therefore,  variations  in 
atmospheric  pressure  occurring  during  the  experiments 
and  or  during  calibrations,  will  influence  the  water 
surface  elevations  measured  by  the  pressure  sensor. 

As  an  excimple,  a variation  of  .5  inch  of  mercury  in 
the  atmospheric  pressure,  will  give  a difference  of 
approximately  17  cm  in  the  mean  water  depth  computed 
with  the  pressure  sensor  data. 

The  difference  between  the  mean  water  level  measured  using 
the  wave  staff  and  pressure  sensor  ranged  between  8 to  20  cm, 
as  given  in  Table  II.  The  approximate  decrease  in  pressure 
head  due  to  the  Bernouli  term  is  given  by  the  variance  of  the 
on-offshore  horizontal  velocity.  The  Bernoulli  contribution 
varied  between  5 to  27  cm,  which  is  a substantial  contribution 
The  magnitude  of  the  Bernoulli  term  varies  with  location  rela- 
tive to  the  breaker  point  and  the  type  of  the  breaker.  Hence, 
it  is  concluded  that  substantial  errors  can  be  introduced  in 
the  mean  water  elevation  in  very  shallow  water  using  pressure 
sensors  due  to  the  Bernoulli  term. 


300  005  010  01 S oao  025  030 

Figure  3.  Mean  water  depth  given  by  the  wave 

gage  and  by  the  pressure  sensor,  19  July 
(evening) 
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B.  SURFACE  ELEVATION  SPECTRUM  CALCULATED  FROM  PRESSURE 

SPECTRUM  USING  LINEAR  TRANSFER  FUNCTION  (STANDARD  TECHNIQUE) 

The  pressure  spectrxim  was  converted  to  a theoretical  wave 
spectrum  for  comparison  with  the  measured  wave  staff  spectrum 
using  the  transfer  function  given  by  linear  theory.  (Equation 
(3)  in  section  IV.)  The  power,  coherence  squared  (henceforth 
referred  to  as  coherence) , and  phase  spectra  of  the  wave  staff 
and  pressure  sensor  are  shown  in  Fig.  4 for  the  19  July, 
evening.  During  the  19  July,  evening  measurements,  the  tide 
was  high  resulting  in  the  wave  sensors  being  located  just  out 
of  the  surf  zone;  hence,  the  waves  are  indicative  of  very 
shallow,  non-breaking  waves. 

The  power  spectra  shows  a narrow-banded  peak  near  a fre- 
quency of  .0938  Hz  corresponding  to  a period  of  10.7  sec.  Also 
evident  are  peaks  at  approximately  .1875  and  .2734  Hz,  which 
appear  to  be  harmonics  of  the  primary  frequency  at  .0938  Hz. 

In  all  the  data_  there  is  a relative  peak  of  energy  at  a 
frequency  of  approximately  .02  Hz,  which  we  attribute  to  surf 
beat.  Above  approximately  .4  Hz,  the  linear  transfer  function 
over-compensates  and  the  surface  elevation  spectrum  calculated 
from  the  pressure  records  diverges  rapidly. 

The  coherence  between  surface  elevation  and  pressure  of 
all  experiments  were  high,  ranging  from  .95  to  .99  in  the  max- 
imum energy  portion  of  the  wave  spectra.  These  extremely  high 
coherences  begin  rapidly  decreasing  at  approximately  .45  Hz. 

The  8 March  coherence  spectrum  is  similar  with  the  exception 
of  the  region  below  .05  Hz.  In  this  region,  the  8 March  data 
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shows  coherence  values  of  .80  to  .85  with  a minimum  of  .68  at 
approximately  .107  Hz. 

The  phase  difference  between  the  pressure  sensor  and  wave 
staff  over  the  highly  coherent  band  of  prominent  wave  energy 
is  always  small  in  the  range  of  0 to  7 degrees.  Because  the 
pressure  sensor  was  located  approximately  40  cm  shoreward  of 
the  wave  staff,  the  wave  staff  phase  always  leads  the  pressure 
sensor  phase.  During  the  8 March  experiment,  the  pressure 
sensor  was  directly  beneath  the  wave  staff  and  the  signals  are 
in  phase.  The  phase  angle  becomes  random  for  the  non-coherent 
region  of  the  spectrum. 

C.  SURFACE  ELEVATION  SPECTRUM  CALCULATED  FROM  PRESSURE 

SPECTRUM  INCLUDING  THE  BERNOULLI  TERM 

As  described  in  section  III,  two  techniques  are  proposed 

2 2 

here  to  include  the  usually  neglected  Bernoulli  term,  u + v , 
in  converting  pressure  to  surface  elevations.  The  spectral 
densities  of  surface  elevations  at  selected  frequencies  for 
the  19  July  evening  data,  calculated  using  the  three  pressure 
correction  techniques  compared  in  Table  III. 

1.  Bernoulli  Term  Calculated  Using  Measured  Velocities 

The  first  technique  employs  the  Bernoulli  term  calcu- 
lated using  the  measured  velocities.  The  first  technique  is 
compared  to  the  standard  technique  in  Table  III;  the  spectra 
are  shown  in  Fig.  5.  In  frequencies  ranging  from  .0156  to 
.0313  Hz  (surf  beat  region),  the  first  technique  has  decreased 


TABLE  III.  SURFACE  ELEVATION  SPECTRA  FOR  19  JULY  EVENING 


(Hz)  at  spectral  peaks;  S = spectral  densities;  n = wave  gage; 
sensor;  1 = first  technique;  2 = second  technique 


F»'i"  rr — ^ 

the  spectral  densities  of  the  converted  pressure  spectrum  to 
values  closer  to,  but  less  than,  the  wave  staff  spectral  den- 
sities. Referring  to  the  spectral  density  ratios  listed  in 
columns  3-5  of  Table  III,  it  is  noted  that  for  the  same  fre- 
quencies, application  of  technique  1 improves  the  spectral 
density  ratios  from  1.036  to  .987,  1.022  to  .980  and  from 

I 

1.074  to  1.021  respectively;  a relative  improvement  in  all 
cases. 

Without  the  introduction  of  the  non-linear  Bernoulli  term, 

the  pressure  spectral  densities  are  consistently  smaller  than 

the  spectral  densities  of  the  wave  gage  spectrum  in  the  band 

of  frequencies  from  .0781  to  .1172  Hz  where  most  of  the  energy 

2 2 

IS  present.  After  the  application  of  the  u + v term  as 
evaluated  by  technique  1,  all  values  show  improvement  with  the 
exception  of  the  frequency  . 1016  Hz  where  the 
ratio  of  .985  was  changed  to  1.017.  The  same  is  generally  true 
for  the  harmonic  spectral  peaks,  with  the  exception  of  the  fre- 
quency .1953  Hz.,  where  the  ratio  of  .977,  the  application  of 
the  first  technique  gives  a ratio  of  1.024  (see  Fig.  6). 

For  the  8 March  data  (Table  IV) , in  the  frequency  band  that 
contains  most  of  the  energy,  the  standard  pressure  spectrum 
spectral  densities  are  uniformly  smaller  than  those  of  the 
wave  staff  spectrum,  with  the  exception  of  the  surf  beat  fre- 
quency of  .0098  Hz.  Application  of  the  first  technique  in  the 
surf  beat  frequency  of  .0098  Hz  increased  the  spectral  density 
ratio  whereas  in  the  first  harmonic  (.1270  Hz)  we  reduced  the 
ratio  .941  to  .854. 


In  the  itiaximuin  energy  band  range,  technique  1 gives  better 
agreement  for  two  of  the  frequencies  (.0586  and  .0781  Hz),  but 
gives  slightly  poorer  agreement  in  the  other  two  frequencies 
(.0684  and  .0879  Hz). 

2 . Bernoulli  Term  Calculated  by  Convolving  Pressure 

Records 

In  the  second  technique,  the  velocity  is  calculated  by 
convolving  the  pressure  record,  using  the  transfer  function 
derived  from  linear  theory  as  given  by  Eq.  12.  A typical 
analogue  record  of  both  measured  and  computed  vector  magni- 
tudes is  shown  in  Fig.  6.  The  velocity  vectors  of  the  meas- 
ured velocity  and  the  computed  velocity  determined  by  convolu- 
tion of  pressure  are  compared  in  Fig.  7;  the  power  spectra  are 
very  similar  and  high  coherence  and  phase  agreement  in  the  band 
of  frequencies  smaller  than  approximately  .3  Hz  are  to  be  noted. 
Hence,  it  is  concluded  that  the  use  of  velocities  calculated 
from  the  convolved  pressure  records  should  give  reasonable 
results  compared  to  direct  measurements. 

For  the  19  July  evening  experiment,  the  improvement 
introduced  by  the  second  technique  is,  in  fact,  always  better 
than  the  improvement  obtained  with  the  first  technique,  with 
the  exception  of  the  frequency  of  .0188  Hz,  (Fig.  8).  For  the 
8 March  data,  compared  with  the  first  technique,  the  improve- 
ment obtained  with  the  second  technique  is  better  in  all  fre- 
quencies with  the  exception  of  the  frequencies  of  .0586  and 
.0781  Hz. 
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3 . Comparing  Variances 


Variances  (and  their  ratios)  calculated  by  integrating 
the  spectra  from  .007  to  .30  Hz,  are  given  in  Table  V.  For  all 
the  experiments  the  variance  of  the  surface  elevation  spectra 
computed  from  the  pressure  record  using  the  conventional  trans- 
fer functic  's  smaller  than  the  variance  of  the  wave  staff 
spectra,  w . the  exception  for  the  20  July  evening  data  where 
it  is  slightly  larger. 

For  the  19  July  evening  data  the  improvement  on  the 
error  was  3.4  and  3.1%  using  first  and  second  technique  respec- 
tively. The  best  improvement  obtained  was  for  the  19  July 
morning  experiment  where  the  improvement  was  7.8%  for  technique 
1 and  9.7%  using  the  second  technique. 

For  all  other  experiments,  the  inclusion  of  the 
Bernoulli  term  resulted  in  an  over  compensation  of  surface 
elevation  variance. 


) 

i 


I 

! 


) 


i 

( 


I 


TABLE  V.  COMPARING  SURFACE  ELEVATION  VARIANCES 


VI.  CONCLUSIONS 
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Field  measurements  of  waves  in  very  shallow  water  just 
outside  the  surf  zone  and  at  the  breaking  were  measured  with 
both  a wave  staff  and  pressure  sensor.  The  surface  elevation 
was  inferred  from  the  pressure  using  the  conventional  linear 
wave  theory  transfer  function  and  a method  which  includes  the 
non-linear  quadratic  velocity  term  in  the  Bernoulli  equation. 
The  velocities  were  measured  using  a flowmeter  (technique  1) 
and  calculated  by  convolving  the  pressure  record  (technique  2) . 
The  surface  elevation  inferred  from  the  pressure  sensor  are 
compared  with  the  wave  staff  records  as  a standard. 

Surface  elevations  inferred  from  the  pressure  sensor  were 
generally  underestimated,  mainly  in  the  region  of  the  crest, 
compared  with  the  same  surface  elevation  measured  by  the  wave 
staff.  Pressure  records  are  more  smoothed  than  wave  staff 
records  and  for  energetic  computational  purposes  the  energy 
of  the  wave  measured  with  the  pressure  transducer  is  generally 
smaller  than  that  measured  with  the  wave  staff. 

The  mean  water  levels  determined  from  the  hydrostatic  pres- 
sure head  were  consistently  less  as  compared  with  the  wave 
staff.  The  difference  can  be  due  to  variations  of  the  baro- 
metric pressure  (which  was  not  accounted  for)  and  due  to  the 
dynaunic  pressure  head  which  acts  to  always  underestimate  the 
mean  water  level.  Direct  estimates  show  that  the  errors  in 
the  mean  water  level  due  to  the  dynamic  pressure  head  ranged 
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from  5 to  13  percent.  Therefore,  the  use  of  pressure  sensors 
to  infer  the  mean  water  level  must  be  interpreted  with  caution. 
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It  is  advisable  to  bury  the  pressure  transducer  in  the 
sea  bed  to  insure  that  it  is  out  of  the  flow  field  and  the 
dynamic  pressure  is  essentially  zero. 

In  the  band  of  frequencies  less  than  approximately  .4  Hz 
the  pressure  sensor  and  wave  gage  spectra  are  highly  coherent 
and  almost  perfectly  in  -phase,  whereas  at  higher  frequencies 
the  spectra  become  less  coherent  and  the  phases  random,  as 
noted  previously. 

Both  proposed  techniques  for  calculating  the  surface  ele- 
vation spectrum  from  the  pressure  records  show  a general  im- 
provement as  compared  with  the  standard  technique  using  a 
linear  transfer  function.  For  technique  1,  improvement  (on 
the  order  of  3.4  to  7.8%  of  the  total  variance  of  the  wave 
gage  spectrum)  was  obtained.  Whereas,  for  the  second  technique 
the  improvement  was  on  the  order  of  3.1  to  9.7%.  The  improve- 
ment of  both  techniques  is  approximately  the  same  with  the 
disadvantage  of  the  first  technique  that  the  velocities  are 
measured  with  flowmeters.  In  technique  2 the  velocities  are 
determined  by  convolving  the  pressure  records  which  is  a great 
simplification  in  an  area  where  the  emplacement  of  sensors  is 
difficult,  such  as  the  surf  zone. 

For  both  techniques,  the  reason  why  the  improvement  is  not 
so  large  can  be  explained  by  three  factors: 

1.  The  irrotationality  assiimption  applied  during  the 
development  of  the  Bernoulli  equation  leads  to 
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problems  when  the  equation  is  applied  to  a rota- 
tional flow.  Breaking  waves  are  characterized  by 
both  rotational  and  highly  non-linear  flow. 

2.  In  cases  where  foam-bearing  waves  are  being  meas- 
ured, wave  heights  measured  by  the  wave  staff  are 
always  in  excess  of  those  indicated  by  the  pressure 
transducer.  The  wave  gage  measures  the  wave  ampli- 
tude as  being  to  the  top  of  the  superimposed  foam, 
whereas  the  added  pressure  at  the  pressure  trans- 
ducer due  to  the  focun  is  relatively  small. 

3.  The  Bernoulli  term  is  only  a second-order  correction. 
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APPENDIX  B 


Date 


Instrument 


Gala* 

(meters) 


Calx** 

(meters/volt) 


8 March 


19  July*** 


20  July**’ 


Wave  Gage 

.707 

.683865 

Press.  Sensor 

-1.404 

2.791088 

Flowm.  X Direct. 

- .009 

3.559351 

Flovmi.  Y Direct. 

- .008 

3.532642 

Wave  Gage 

3.455 

- .00122 

Press.  Sensor 

-1.040 

‘ . 006841 

Flowm.  X Direct. 

0.0 

.002251 

Flowm.  Y Direct. 

0.0 

.002261 

Wave  Gage 

3.530 

- .001235 

Press.  Sensor 

-1.040 

.006841 

Flowm.  X Direct. 

o 

• 

o 

. 002251 

Flowm.  Y Direct. 

0.0 

.002261 

Calibration  Additive  Factor 
Calibration  Multiplicative  Factor 

Due  computer  record  technique  the  true  Calx  of  the 
instruments  have  been  multiplied  by  a factor  of  5/2048, 
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SPECTRAL  DENSITIY!hx2-5EC 


POWER,  COHERENCE  AND  PHASE  SPECTRA  FOR  20  JULY  MORNING 


STANDARD  TECHNIQUE 


POWER,  COHERENCE  AND  PHASE  SPECTRA  FOR  20  JULY  EVENING 


TABLE  VI.  SURFACE  ELEVATION  SPECTRA  FOR  19  JULY  MORNING 
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(Hz)  at  spectral  peaks;  S = spectral  densities;  n = wave  gage 
sensor;  1 = first  technique;  2 = second  technique 


TABLE  VII.  SURFACE  ELEVATION  SPECTRA  FOR  20  JULY  MORNING 


(Hz)  at  spectra  peaks;  S = spectral  densities;  n = wave  gage 
sensor;  1 = first  technique;  2 = second  technique 


TABLE  VIII.  SURFACE  ELEVATION  SPECTRA  FOR  20  JULY  EVENING 


(Hz)  at  spectral  peaks;  S = spectral  densities;  n = wave  gage 
sensor;  1 = first  technique;  2 = second  technique 
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